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Abstract. We revisit the classical problem of the buckling of a long thin axially compressed cylindrical shell. 
By examining the energy landscape of the perfect cylinder we deduce an estimate of the sensitivity of the shell to 
imperfections. Key to obtaining this is the existence of a mountain pass point for the system. We prove the existence 
on bounded domains of such solutions for all most all loads and then numerically compute example mountain pass 
solutions. Numerically the mountain pass solution with lowest energy has the form of a single dimple. We interpret 
these results and validate the lower bound against some experimental results available in the literature. 
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1. Introduction. 

1.1. Buckling of cylinders under axial loading. A classical problem in structural en- 
gineering is the prediction of the load-carrying capacity of an axially- loaded cylinder. As well 
as being a commonly used structural element, the axially- loaded cylinder is also the archetype 
of unstable, imperfection-sensitive buckling, and this has led to a large body of theoretical and 
experimental research. 

In the decades before and after the second world war a central problem was to understand 
the large discrepancy between theoretical predictions and experimental observations, as shown in 
Figure 1.1. A variety of different explanations has been put forward, but with the experimental 
work of Tennyson [30] and the theoretical work of Almroth [1] it became clear that this discrepancy 
is mostly due to imperfections in loading conditions and in the shape of the specimens. Further 
experimental and theoretical work by many others has confirmed this conclusion [14, 33, 36]. 

For near-perfect cylinders the linear and weakly non-linear theories (see Section 1.2) ade- 
quately describe the experimental buckling load^ and the deformation just before failure (see 
e.g. [4]). Cylinders used in practical applications, however, are far too imperfect for this weakly 
non-linear theory to apply, and from a practical point of view the problem of predicting the failure 
load is still open. 

In fact there is good reason to believe that it will never be possible to accurately predict 
failure loads for cylinders that are used in practice. For simple materials, such as metals, it is 
believed that current numerical methods can describe the local material behaviour with enough 
accuracy that correct prediction of the complete behaviour of the cylinder — including its failure — 
is feasible. This could be achieved provided the geometrical and material imperfections as well as 
the loading conditions are determined in sufficient detail. The difficulty lies in the qualifier "in 
sufficient detail", since an extremely accurate measurement of geometric imperfections would be 
necessary [4] , and in the design phase both the loading conditions and the geometrical and material 
imperfections in the finalized product are only known in vague terms. Therefore, in recent decades 
the attention of theoretical research has turned to characterizing the failure load in weaker ways, 
preferably in the form of a lower (safe) bound. 

1.2. Characterizing sensitivity to imperfections. Viewed as a bifurcation problem, the 
buckling of the cylinder is a subcritical symmetry-breaking pitchfork bifurcation (Figure 1.2). 
Generically, imperfections in the structure eliminate the bifurcation and round off the branch of 
solutions^, resulting in a turning-point at a load Pimp strictly below the critical (bifurcation) load 
Per of the perfect structure. In an experiment in which the load is slowly increased, the system 
will fail {i.e. make a large jump in state space) at load Pimp- 
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Fig. 1.1. Experimental data from various research groups, all representing failure loads of axially -loaded 
cylinders. The horizontal axis is the ratio of cylinder length and wall thickness; the vertical axis is the ratio of the 
failure load and the theoretical critical load as predicted for perfect cylinders. Note that all tested cylinders fail at 
loads significantly lower than that predicted by theory; in some cases failure occurred at less than one-fifth of the 
theoretical load-carrying capacity. Power-law fitting lines are added to emphasize the dependence of failure load on 
geometry. The data are from [10, 7, 5]. 
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Fig. 1.2. Illustration of perfect and imperfect bifurcation curves. 

Again, if the imperfections in geometry and loading are fully known, then calculation of 
Pimp is a practical problem rather than a theoretical one, and we do not address this problem 
here. For the more difficult question of characterizing Pimp under incomplete information various 
strategies have been proposed. A classical line of thought originates with Koiter [22], in which 
the imperfections are chosen a priori within certain finite-dimensional sets. Common choices are 
the sets spanned by an eigenvector at bifurcation, or by the eigenvectors associated with the first 
n bifurcations. This approach might be termed weakly non-linear, as it is based on an expansion 
of the energy close to the bifurcation point, in the directions suggested by the bifurcation point 
itself. It gives predictions that are correct if the imperfections are very small — much smaller than 
those encountered in practice. 

Since the a priori choice of imperfections is a weak point of this method, a natural step is 
to optimize over all possible perturbations. Demi and Wunderlich pioneered this approach, in 
which a numerical algorithm is used to find a "worst geometric imperfection" [9]. This "worst 
imperfection" is defined as the one that produces a turning-point of minimal load. Some constraint 
on the magnitude of allowable perturbations is necessary, of course, to prevent running a steam 
roller over the cylinder being interpreted as an admissible imperfection. The authors of [9] and [35] 
first suggest constraining the P^-norm of the perturbation displacement, but they immediately 
replace the P'^-norm by an P^-norm for computational convenience. 

This method has an interesting aspect that is often glossed over in the engineering literature. 
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By definition ttie failure load obtained by this method is a lower bound for the failure load of all 
systems that have perturbations of less or equal magnitude. The measure of magnitude is defined 
by the choice of constraint. Therefore the choice of constraint on the imperfections is a critical 
one, since it implicitly defines a class of imperfections that produce the same or higher failure 
load. 

1.3. Main results. In this paper we follow a related, but distinct, line of reasoning. Instead 
of studying actual behaviour of imperfect cylinders, we deduce an estimate of the sensitivity to 
imperfections from the energy landscape of the perfect cylinder. The final result is a lower bound 
on the failure load similar to above, and the approach gives additional insight into the problem. 

The key result is the existence of a mountain-pass pointy an equilibrium state that sits astraddle 
in the energy landscape between two valleys; one valley surrounds the unbuckled state, and the 
other contains many buckled, large-deformation states. 

This mountain-pass point has a number of interesting properties: 

1. It has the appearance of a single- dimple solution^ a small buckle in the form of a single 
dent (see Figure 1.3 (a)). Single-dimple deformations have appeared in the engineering 
literature in a number of different ways (see Section 6), but theoretical understanding of 
this phenomenon is still lacking. Localization (concentration) of deformation is known 
commonly to appear in extended structures [21], and in the cylinder localization in the 
axial direction has been studied theoretically and numerically [20, 24, 25]. Whether local- 
ization is possible in the tangential direction has been an open problem for some time; it is 
interesting that our simulations for the perfect structure show solutions that are localized 
in both axial and tangential directions. 

2. Like all mountain-pass points this single-dimple solution is unstable, in the sense that there 
are directions in state space in which the energy decreases. In one direction the dimple 
roughly shrinks and disappears, and in the opposite direction it grows and multiplies 
(Figure 1.3 (b-c)). It is remarkable, however, that our numerical results indicate that the 
single-dimple solution has an alternative characterization as a constrained global minimizer 
(a global minimizer of the strain energy under prescribed end shortening). 

3. The equations can be rescaled so that the only remaining parameters are the load level 
and the domain. The geometry of the mountain-pass solution we calculate even appears 
to be independent of the domain size. 

This mountain-pass point is central in an estimate of the sensitivity to perturbations. For 
the system to escape from the neighbourhood of the unbuckled state, it must possess at least the 
energy associated with this mountain-pass point. The mountain-pass energy level is therefore an 
indication of the degree of stability of the unbuckled state. Implicitly it defines a class of perturba- 
tions for which the unbuckled state is stable. This approach is related to the "perturbation-energy" 
approach first suggested by Kroplin and co-workers [23, 11, 32], but differs in some essential points; 
see the discussion in Section 5. 

1.4. Methods. We use both analytical and numerical methods. In Section 2 we introduce 
the Von Karman-Donnell equations, which form the basis of this paper, and rescale them in an 
appropriate manner. In Section 3 we present the functional setting that we use, show that the 
energy functional has the geometry associated with a mountain pass, and prove the existence 
of mountain-pass points (Lemma 3.5). There are certain interesting technical issues. By their 
localized nature, single-dimple solutions are most naturally defined on an unbounded domain; 
however, we are only able to prove existence of mountain-pass points on bounded domains, and 
consequently we work on finite domains that become large in the limit of thin shells. Similarly, 
the non-coercive nature of the energy functional implies that we prove existence of mountain-pass 
points for almost all load levels (see Lemma 3.4). 

In Section 4 we turn to numerical investigation. We use a variety of different algorithms 
to find solutions of the discretized Von Karman-Donnell equations. With a discrete mountain 
pass algorithm we find solutions that are, by construction, mountain-pass points. The solution 
of Figure 1.3 (a) was found in this manner. With a constrained gradient fiow we also find local 
minima of the strain energy under prescribed end shortening. Some of these solutions appear 
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Fig. 1.3. Part (a) shows a numerical computation of a solution w-^p that is a mountain-pass point of the 
energy Fx for a load A = 1.5. We show both the graph of w-^p(x,y) as well as its rendering on a cylinder. At 
wmp there exist two directions in the state space X in which the energy Fx decreases. By perturbing w-^p in these 
directions and following a gradient flow of Fx we move away from w-^p. In one direction the dimple shrinks to 
zero (not shown) and in the opposite it grows in amplitude and extent (figures (b) and (c)). 

to coincide with those found by the discrete mountain pass algorithm, and the mountain pass 
solutions are stable under this gradient flow. These observations lead us to conjecture that the 
global mountain pass solution is also a global constrained minimizer of the strain energy. By a 
constrained version of the discrete mountain pass algorithm we also find critical points of higher 
Morse index. 

Section 5 is devoted to an interpretation of these results in the context of imperfection sensi- 
tivity, as mentioned above, and in Section 6 we wrap up with the main conclusions. 

2. The Von Karman-Donnell equations. We consider a cylindrical shell of radius i?, 
thickness Young's modulus and Poisson's ratio that is subject to a compressive axial 
force P. In dimensionless form the Von Karman-Donnell equations are given by 

e^/\^w + AUJ^^ - ~4>^^ - 2[w,'4>] = 0, (2.1) 
A2^+[1U,1U]+1U^^ = 0, (2.2) 

where subscripts x and y denote differentiation with respect to the spatial variables, and the 
bracket [•, •] is defined as 

[(2, h\ = ~CLxxbyy H~ "^^yy^xx ^xy^xy 

The function w is the inward radial displacement measured from a unbuckled (fundamental) state, 
(j) is the Airy stress function, = t^(1927r^i?^(l — z/^))~^ and the nondimensional load parameter 
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is given by A = P{S7r^ERt)~^ . The unknowns w and (j) are defined on the two-dimensional spatial 
domain {—i^i) x (—1/2, 1/2), where x G {—i^i) is the axial and y G (—1/2, 1/2) is the tangential 
coordinate. Since the ^/-domain (—1/2,1/2) represents the circumference of the cylinder, the 
functions w and ^ are periodic in y; at the axial ends x G {—i^i} they satisfy the boundary 
conditions 

Wx = {Aw)x = (l^x = {^^)x = 0. 
In Appendix B we discuss the derivation of this formulation. 

Equations (2.1-2.2) can be rescaled to be parameter-independent in all but the domain of 
definition and load parameter, by substituting 

w ^ ew^ (j) ^ £^(j)^ X ^ e^l^x^ y ^ s^^^y-, A ^ eX. (2.3) 

so that the equations become 

A^^ + Xwxx -(l)xx-2 [w, (j)] = 0, (2.4) 

/^^(j)^[w,w\ ^-Wxx = 0. (2.5) 

The domain of definition of w and (j) is now 

n := (-fe-l/^££-l/2) X (-|£-^/^ ^^"'/'), 

which expands to as 6 ^ 0. When not indicated otherwise we choose the aspect ratio 2-^ = 1; 
in Section 4.3 we comment on the infiuence of domain size and aspect ratio. 
The boundary conditions for w and (j) now are 

w is periodic in y^ and Wx = {Aw)x = at x = i^e""*^/^, (2.6a) 
(j) is periodic in y^ and (j)x = {^(l))x = at x = i^e""*^/^. (2.6b) 

Equations (2.4-2.5) are related to the stored energy E and the average axial shortening S 
given by 

E{w) -.= 1 [ {Aw^ + A4>^) , and S{w) := ^ / wl (2.7) 

Note that the function (j) in (2.7) is determined from w by solving (2.5) with boundary condi- 
tions (2.6b). Furthermore, solutions of (2.4-2.5) are stationary points of the total potential 

Fx{w) :=E{w)-XS{w) (2.8) 

and are also stationary points of E under the constraint of constant S. We use both properties 
below. 

3. The mountain pass: overview. We briefiy recall the general context of the Mountain- 
Pass Theorem of Ambrosetti and Rabinowitz [2]. Let / be a functional defined on a Banach space 
X, and let wi, W2 be two distinct points in X. Consider the family F of all paths in X connecting 
Wl and W2 and define 

c= inf max/(K;) , (3.1) 

that is the infimum of the maxima of the functional / along paths in F. If c > max{/(K;i), I{w2)}^ 
then the paths have to cross a "mountain range" and one may conjecture that there exists a critical 
point wmp of / at the level c, called a mountain pass point. 

We will apply this idea to the Von Karman-Donnell-equations in the following way. We take 
for / the total potential Fx (see (2.8)) at some fixed value of A, and for the end point Wi the 
origin. We will obtain a mountain-pass solution by the following steps: 
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MPl. We first show that k;i = is a local minimizer: there are g^a > such that Fx{w) > a 

for all w with \\w\\x = Q (Lemma 3.1); 
MP2. If £ is small enough, then there exists W2 with Fx{w2) < (Corollary 3.3). 
MP3. Given a sequence of paths 7^ that approximates the infimum in (3.1), we extract a (Palais- 
Smale) sequence of points Wn ^ Jn^ each one close to the maximum along 7^, and show 
that this sequence converges in an appropriate manner (Lemmas 3.4 and 3.5). 
In this way it follows that there exists a mountain-pass critical point w with Fx{w) = c, 
provided that e is sufficiently small (or that the domain is sufficiently large). For technical reasons 
(lack of coerciveness of the functional Fx) this procedure can be performed only for almost all 
< A < 2 (see Lemma 3.4). 

In the rest of this section we detail the steps outlined above. 

3.1. Choice of spatial domain. We are interested in mountain-pass solutions of the system 
of equations (2.4-2.5) that are quasi-independent of the domain size in the sense that they 
converge to a non-trivial solution on as £ ^ 0. This point of view suggests considering the 
problem on the whole of M? rather than on a sequence of domains of increasing size; however, 
there are two reasons not to do this. To start with, the numerical calculations described below are 
necessarily done on a bounded domain; more importantly, for the proof of existence of mountain- 
pass points boundedness of the domain is necessary. For these reasons we concentrate on bounded 
domains, while keeping the context of the unbounded domain in mind. 

3.2. Functional setting and linearization. We introduce a functional setting for the 
functions w that is suggested by the linearization of the stored energy functional E. Writing 
^ = 01 + 02, where 

A'^(l)i = -Wxx and A^02 = -[^,^], (3.2) 
we can expand the energy functional E as 

Eiw) = i / Au;2 + 1 /■ + /" A0iA</.2 + i / Acj^l (3.3) 

Since 02 is quadratic in the second derivative of E is given by 

d^E{0)-u-v= [ AuAv^ [ A0^A0^, 

where 0^^'^ are obtained from u and v by (3.2). Inspired by this linearization of E we define 

X = G H'^{n) : V^^(±^£"^/^ •) = 0. is periodic in y, and J = o| 

with norm 

\\w\\l= [ ( Aw'+A^l), 

JQ 

where 0i G H'^{Q) is the unique solution of 

A^0i = —Wxx-, 01 satisfies (2.6b), and / 0i = 0. 

Jq 

This norm is equivalent to the i^^-norm on the set X, and with the appropriate inner product the 
space X is a Hilbert space. 

We now address the requirements of the mountain pass theorem mentioned above in MPl- 
MP3. 
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3.3. The origin is a local minimizer. The norm in X is related in a natural manner with 
the shortening 5, as is demonstrated by the (sharp) estimate 



/ wl = - f ww^^ = f wA^(l)i = f AwA(l)i <l f Aw^ f A(l)l = l \\w\ 
JQ Jn Jn Jn ^ Jq ^ Jq ^ 



2 

X ' 



(3.4) 

This inequality strongly suggests that for A < 2 the origin is a strict local minimum for the 
functional Fx{w) = E{w) - XS{w). 
Lemma 3.1. 

1. There exists a constant C , dependent on Q., such that 



<C\\w\\%. 



(3.5) 



2. For any < A < 2 there exists ^ > such that 

■mi{Fx{w) : ||u;||^ = ^,} >0. 

Proof. Split = 01 + 02 as in (3.2), and note that the function V(/)i is bounded in by the 
Sobolev imbedding {il; e : / = 0} L"^: 

The third term on the right-hand side of (3.3) can now be rewritten as 



A(/)iA02 
which we estimate by 



W,W\ = / (l)l{w^xWyy - Wl ) = 



WxWxy - (l^lxWxWyy), 



A0iA02 



< 2 ||V0i||^. 11^.11^. IIA^II^. < C \\wxh, \\Aw\\l, < CVSHMx . 



thus proving the first part of the Lemma. 

Since A < 2, choose < £» < (2 - A)/2C and define t] = ^{l-Cg- A/2) > 0. Then on the set 
C = {w : \\w\\^ = g}, using (3.4), we find that 

Fxiw) = E{w) - XSiw) 



>-\\w\\i,-C./SMMx-^S{w) 
-2 ^2 4 



□ 

Part 2 of this lemma implies that by choosing wi to be the origin we have shown condition 
MPl. 

Remark 3.1. Although the inequality (3.4) suggests that the origin should be a local minimizer 
for any domain 1], bounded or not, the proof above only applies to bounded domains. F. Otto 
has constructed a proof of this result that is valid on any domain (personal communication). 
Interestingly, this proof uses not only the cubic energy term A0iA02 but also the quartic term 
A(/)2, and appears to break down without this latter term. 
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3.4. Periodic solutions exist with negative Fx* To satisfy MP2 we show in this section 
that for any A > functions w G X exist for which Fx{w) = E{w) — XS{w) < 0. To do this we 
construct a sequence of functions ws with specific scahng properties: 

Lemma 3.2. There exists a sequence of functions ws, 1-periodic on such that 

I ^1, I Awl = 0{S-'), 

[-l/2,l/2]2 [-l/2,l/2]2 

and J = 0((52-^) as 5^0, 

[-l/2,l/2]2 

for any a > 0. Here the functions ws and(j)s solve equation (2.5) with periodic boundary conditions. 
In addition, ws and (j)s satisfy the boundary conditions (2.6) on the boundary o/ [—1/2, 1/2]^. 

The proof, given in the Appendix, is inspired by the so-cahed Yoshimura pattern [37], a 
folding pattern by which a fiat sheet of paper, or a cyhndrical sheet of thin material, can be folded 
into a macroscopically cylindrical structure with zero Gaussian curvature but locally infinite total 
curvature (Figure 3.1). The functions ws are smoothed versions of the Yoshimura pattern, adapted 
to the geometrically linear setting of the von Karman-Donnell equations, and S measures the width 
of the fold. 

Corollary 3.3. 

1. Fix A > 0. If e is sufficiently small, then there exists w e X such that Fx{w) < 0. 

2. Fix e sufficiently small. Then there exists Xo{£) G [0,2) such that for all A > Aq there 
exists w G X with Fx{w) < 0. 

Proof By scaling the functions ws of Lemma 3.2 the claims can be fulfilled: let S = and 

set 

Ws{x,y) = £"^^^2/3(x£:^/^?/£:^/^), 4>s{x,y) = £"^^^2/3 (xe^/^, ^/e^/^); 
then G X, and (2.5) is invariant under this scaling; in addition, choosing a = 1/6, 




Therefore lim^^o = 0, proving the first claim. For the second claim we fix e such that < 2; 
then for ah A > Q^, Fx{ws) < 0. □ 




Fig. 3.1. Yoshimura folding pattern. 
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3.5. Convergence of selected sequences. For given A G (0, 2) and for sufficiently small 
£ > 0, the two previous sections provide two points: the origin wi = that satisfies MPl and a 
point W2 with Fx{w2) < 0, such that 

c(A) := inf maxFA(^^;) > 0, (3.6) 

where V is the set of curves connecting and 

r = {7 G C([0, 1]; X) : 7(0) = 0, 7(1) = W2}, 

(actually, V depends on A through the dependence on but W2 can be taken independent of A 
in a neighbourhood of a given A G (0, 2)). 

We were unable to prove the classical Palais-Smale condition, which reads 

For any sequence Wn ^ X such that Fx{wn) c and F^{wn) ^ in X', there 
exists a subsequence that converges in X. 
The difficulty lies in the lack of coerciveness of the functional Fx: the quotient Fx{w)/ \\w\\\ is 
not bounded away from zero, implying that Palais-Smale sequences may be unbounded in X. 

The "Struwe monotonicity trick" [29] provides a way of proving the boundedness of Palais- 
Smale sequences for at least almost all A G (0,2). The pertinent observation is that for fixed 
F\{w) is decreasing in A; by consequence c(A) is a decreasing function of A and therefore differen- 
tiable in almost all A G (0, 2). If ^{t) is the highest point of a near-optimal curve 7 at some Aq, then 
c'{\q) should be close to —S{^{t)). Finiteness of c' at Aq thus implies that near- mountain-pass 
points have bounded and this additional information suffices for the construction of bounded 
sequences: 

Lemma 3.4. Let A G (0, 2) he such that c'{\) exists. Then there exists a hounded Palais-Smale 
sequence Wn, i.e. a sequence that satisfies 

1. Wn is hounded in X; 

2. F'^{wn) -^^ inX' and Fx{wn) ^ c(A); 

3. there exists a sequence of curves (7n) C F such that Wn G 7n([0, 1]) and max^^[o 1] Fx{jn{t)) 
c(A). 

In [28] this same argument was used to study mountain-pass points for the related one- 
dimensional functional 



Jxiu) 



where F is a non-negative double- or single-well potential. The proof of Lemma 3.4 is a word-for- 
word repeat of the proof of [28, Proposition 5], and we omit it here. 

Lemma 3.5. The sequence Wn given hy Lemma 3.4 is compact in X, and a suhsequence 
converges to a stationary point w e X of Fx. 

Strictly speaking the stationary point given by this lemma may not be a mountain-pass point 
itself, in the sense that there may not be a curve 7 G F of which w is the highest point. Property 3 
of Lemma 3.4, however, states that w has an approximate mountain-pass character. 

Proof We extract a subsequence that converges weakly in X and strongly in and L"^ to 
a limit w. Defining and (j)2n by (3.2), we find that the right-hand sides in (3.2) are bounded 
in I/^ and L^, and therefore that (j)in and (j)2n converge strongly (up to extracting a subsequence, 
which we do without changing notation) in to functions 0i^2- Both functions (f)i^2 are again 
related to w by (3.2); for 02 this follows from remarking that for given ( G C^(^l), 

lim / C[Wn,Wn] = lim / Wn[C,Wn] = / w[C,w]= / C[w,w], 

so that the right-hand side converges in the sense of distributions. Similarly it follows from the 
strong ii'^-convergence of (j)ri = </>in + 4^2n that 

lim / (l)n[Wn,W - Wn] = lim / Wn[(l)n,W - Wn] = 0. 
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as 



To show that Wn converges strongly in X, note that the derivative F'^{wn) can be characterized 

F'^{Wn)-V= / /^Wn^V - I (j)n{Vj:j: ^ 2[Wn, v]) - \ I Wnx^x- 

Jn Jn Jq 



We now calculate 



lim < / Aw'^ — / Aw'^ > = lim / AwnA{w — Wn) 

'l(^n) • - ^n) + / (l)n{{w - Wn)xx ^^[Wn,W - Wn]) 

Jn 



lim iF 



= 0. 



The strong convergence of Wn in X now follows from the uniform convexity of X. □ 
4. Numerical Results. 

4.1. Description of the algorithm. Our goal in this section is to find, numerically, critical 
points of Fx. Although we will focus on mountain-pass points described above and sketch the 
method used to find them, numerical approximations of other critical points of Fx will be shown 
as well. More details on all the numerical methods used are given in a companion paper [17]. 

In order to employ the mountain-pass algorithm we discretize (2.4-2.8) using finite differences. 
The algorithm was first proposed in [8] for a second order elliptic problem in ID. It was later used 
in [18] for a fourth-order problem in 2D. 

The main idea of the algorithm is illustrated in Figure 4.1. We take a discretized path 
connecting wi = with a point W2 such that Fx{w2) < 0. After finding the point Zm at which Fx 
is maximal along the path, this point is moved a small distance in the direction of the steepest 
descent of Fx. Thus the path has been deformed and the maximum of Fa lowered. This deforming 
of the path is repeated until the maximum along the path cannot be lowered any more: a critical 
point Wy^p has been reached. 




Fig. 4.1. Deforming the path in the main loop of the mountain pass algorithm: point Zm is moved a small 
distance in the direction —\/Fx{zm) CLnd becomes z^"^ . This step is repeated until the mountain pass point w-^p 
is reached. 

Figure 4.2(a) shows a numerical solution of (2.4-2.5) obtained by the this algorithm with 
A = 1.1. The graph shows the radial displacement as a function of x and y. Rendered on a 
cylinder this solution represents a single dimple as can be seen below the graph. 
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Fig. 4.2. Numerical solutions found using the (constrained) mountain pass algorithm, constrained steepest 
descent method, and the Newton algorithm. The figures show both the graph ofw(x,y), as well as its rendering on 
a cylinder. 

We restrict our computations to functions that are even about the x- and ?/-axes, i.e. to the 
subspace 

y = {iIj eX : y) = y), y) = tlj{x, -y)} , 

thus reducing the computational domain ft to one quarter, e.g. (0, x (0, ^e"^/^). The 
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boundary conditions (2.6a) then become 

W:, (Ai^)^ = for X G {0, and Wy = {/^w)y = for G {0, 

This symmetry assumption has many numerical advantages, but on the other hand it a priori 
excludes solutions that do not belong to y. 

For the mountain-pass algorithm we always use the unbuckled state = as the first 
endpoint of the paths. The choice of the second endpoint W2 has a non- trivial influence on the 
solution to which the mountain-pass algorithm converges. Corollary 3.3 guarantees the existence 
of ^L'2 G ^ with Fx(w2) < 0; in the numerical implementation, however, we found such a k;2 by a 
steepest descent method (rather than taking the function constructed in the proof of Lemma 3.2): 
starting from a function wq that has one peak located in the center of the domain 1], we solved 
the initial value problem 

j^w{t) = -VFx{w{t)), w{0) = wo, (4.1) 

on an interval (0, T) until Fx{w{T)) < 0. We then defined W2 = w{T). 

A different choice of W2 (or, more precisely, of the starting point wq of (4.1)) can lead to 
a different solution of the problem, as Figure 4.2 (b,f) shows. Here wq was chosen to have two 
peaks with centers on the axes x = and ^ = 0, respectively. The algorithm then converged to a 
numerical solution with two dimples in the circumferential and axial direction, respectively. 

Note that the numerical solution Wy^^ selected by the mountain-pass algorithm has the 
mountain-pass property in a certain neighbourhood only: there exists a ball Bp{wy^Y>) ^^o 
points Wi,W2 G Bp{wy^p) such that 

^a(^mp) = max Fa (^) > max{Fx{wi), Fx{w2)} , 

where F is the set of curves in Bp{wy^p) connecting wi and W2. The reason for this is that the 
algorithm deforms a certain initial path connecting wi and W2 which is fixed. In order to recover 
the global character one would need to run the algorithm for all possible initial paths. 

The rest of the numerical solutions shown in Figure 4.2 were obtained under a prescribed 
value of shortening S by the constrained steepest descent method and the constrained mountain 
pass algorithm [16, 17]. 

4.2. Calculation of Fx{w^p). In the preceding sections we have showed that 

1. for a sufficiently large domain 1], a function W2 on Q exists with Fx{w2) < 0; 

2. for each such function W2 and for almost all < A < 2, a mountain-pass solution Wy^p = 
Wy^p{X^Q^W2) exists. 

Different end points W2 may give rise to different mountain-pass points, as we have observed in the 
numerical experiments described above. We therefore define the mountain-pass energy function 
V on (0,2) by 

V{\n) :=mf{Fx{wy,p{\n,W2)) :Fa(^2) <0}. (4.2) 

For a given A, the value of F(A, is the lowest height (or energy level) at which one may pass 
from the neighbourhood of the origin to any point with negative total potential Fx. We now derive 
some of its properties and calculate it numerically. 
Lemma 4.1 (Properties of V(A,1])). 

1. For sufficiently large ft there exists Xo{ft) > such that V(A, 1]) < oo for almost all 
AG (Ao,2); 

2. V is a decreasing function of A; 

3. For sufficiently large there exists c{rt) > such that 

v{x,n) < c{2-xf 



for sufficiently small 2 — A > 0. 
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Proof. Part 1 is a reformulation of the main result of Section 3, making use of Corollary 3.3. 
For part 2 we remark that for each fixed Fx{w) is a decreasing function of A; the infimum of a 
set of decreasing functions is again decreasing. 

For part 3, let us set 

E{w) = E2{w) + E^{w) + E^{w), 

where 

E2{w) := i ||«;||^ = i /" (A«;2 + Esiw) := f ^4>iM2, and Ei{w) :=]- j ^4>l 

where (j)i and 02 are determined from w by (3.2) (see also (3.3)). Note that En has homogeneity 
n, i.e. Eni/iiw) = ii^En(w). 

A classical result in the engineering literature of cyhnder buckling (see e.g. [20]) states that 
there exists a periodic function w on such that 

E2{w) = 2S{w) and E^^iw) + E^{w) < 0. 

Here and below we consider the integrals that define En{w) and S{w) as taken over a single 
period cell. For sufficiently small 2 — A > the inequality above gives that Fx{w) = (2 — X)S{w) + 
Es^w) + E4{w) is negative, implying that w is an admissible end point W2 for the definition (4.2) 
of y(A, 1]), and the connecting hue segment {fiw : < /i < 1} therefore an admissible curve in F. 
Consequently 

V(A) < sup Fxii^w) = sup /i^(2 - X)S{w) + ij^E^{w) + ii^E^(w) 

0</i<l 0</i<l 

The supremum on the right-hand side is obtained at 
_ 3IE3HI [ / 32(2-A)gH.E4H 

implying that the claim holds for periodic functions. The generalization to non-periodic functions 
on large domains Q {i.e. for small s) is made by filling the domain with a large number of periodic 
cells of the function and connecting the function smoothly to the boundary of 1]. □ 

Figure 4.3 shows graphs of the mountain-pass energy V(A, ft) computed for various sizes of 
domain For each domain, the mountain-pass algorithm was employed to compute w-^p for 
several values of A. These mountain-pass solutions were then continued in A using numerical 
continuation. 

4.3. Influence of the domain. The localized nature of the solutions calculated above sug- 
gests that they should be independent of domain size, in the sense that for a sequence of domains 
of increasing size the solutions converge (for instance pointwise on compact subsets). Such a con- 
vergence would also imply convergence of the associated energy levels. Similarly, we would expect 
that the aspect ratio of the domain is of little importance in the limit of large domains. 

We have tested these hypotheses by computing mountain-pass solutions on domains of different 
sizes and aspect ratios. Generally solutions on different domains compare well; the maximal 
difference in the second derivatives of w is two or three orders of magnitude smaller than the 
supremum norm of the same derivative (the details of this comparison are given in [17]). 

Here we only include a calculation of the mountain-pass energy level V{X, Q) for domains Q 
of different aspect ratio and size (Figure 4.3). 

The comparison of solutions computed on different domains and their respective energies 
suggests that for each A we are indeed dealing with a single, localized function defined on R^, of 
which our computed solutions are finite-domain adaptations. In the rest of this paper we adopt 
this point of view, and consequently we will write V{X) instead of ^(A, Q). 

A consequence of this point of view is that dimples in cylinders with different geometric 
parameters are mapped to the same rescaled solution, or equivalently, that the same single-dimple 
solution of (2.4-2.5) corresponds to differently-sized dimples on an actual cylinder, as a function 
of the parameters (Figure 4.4). 



I = ^tItt(2-A)+o(1) asA-2, 

J 3|£:3H| 
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Fig. 4.3. Left: the mountain-pass energy V{X) found numerically for various domains Q. Right: solid line 
shows the same computation as on the left, just plotted for the norm of Wy^-p{X) squared. The dashed curve was 
obtained by continuation of the solid curve; the solutions on the dashed curve do not represent mountain-pass 
points any more. 




Fig. 4.4. The numerical mountain-pass solution w-^-p of the scaled equations (2.4-2.5) for a given value of X 
rendered on two cylinders of the same radius R but a different thickness t: (a) t/R = 0.003; (b) t/R = 0.04. 

5. Interpretation: imperfection-sensitivity. We now turn to the relevance of the moun- 
tain pass in the context of a loading problem. This relevance can best be understood in the context 
of imperfections in the loading conditions (rather than geometric imperfections) such as in the 
case of a (small) lateral loading. 

Under a small lateral load an equilibrium wq^ a local minimum of the functional Fx, may be 
perturbed into an equilibrium wq of a perturbed functional Fx. Since wq is a local minimum, 
^a('^o) > Fxi'^o)^ i'^' with respect to the unperturbed system wq has a higher total potential 
than Wq. The level of Fx that is reached is a measure of the magnitude of the imperfection — a 
different measure than is commonly used, but one that has distinct advantages. 

By definition the number V{\) is the lowest energy level at which it is possible to move between 
the basins of attraction of wi and W2 (Figure 5.1). If the loading imperfection is interpreted, as 
above, as a mechanism capable of maintaining the system at a higher energy level than that of the 
neighbouring fundamental minimizer, then the number V^(A) is critical: as long as the imperfection 
is so small that the energy is never raised by more than V^(A), the new stationary point will be 
part of the same basin of attraction as wi. For larger imperfections, however, it becomes possible 
to leave the fundamental basin of attraction, resulting in a large jump in state space. 

This line of reasoning provides a natural measure of the size of imperfections, namely the 
maximal increase in energy (in the perfect structure) that an imperfection can achieve. It also 
provides a natural measure of the stability of the unbuckled state, since a higher mountain-pass 
energy level implies a larger class of loading imperfections under which the state remains in the 
fundamental basin of attraction. This observation allows us to connect systems with different 
geometrical characteristics, and compare their relative sensitivity to imperfections. 
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5.1. Calibrating the mountain-pass energy. Comparing cylinders of varying geometry 
requires a common measure of imperfection sensitivity It is not a priori clear which measure to 
take; e.g.^ one might consider either the mountain-pass energy itself or the average spatial density 
of this energy which will result in different comparisons for cylinders of different wall volume. 
Here we choose to rescale the mountain-pass energy level by the other energy level present in the 
loaded cylinder: the energy that is stored in homogeneous compression of the unbuckled shell. 

This calculation can be done in two slightly different ways. The first and most straightforward 
is to rescale the dimensional mountain-pass energy (see (B.ll) and (2.3)) 

U.^Em^e^ViX) = 3(3(,_^,))3/.^ ^(A), 

by the elastic strain energy stored in the full length of the compressed cylinder of length L, 

L ^2 ^ ^t'EL ^2 



Ai^ERt 12(1 - 

to give an energy ratio, or a rescaled mountain-pass energy level, 

1 t ViX) 

27rV3(l - L A2 



(5.1) 



From this expression and the calculation shown in Figure 4.3 curves may be drawn in a plot of 
load versus the ratio L/t (see Figure 5.2). Note that to obtain this figure from Figure 4.3 the 
curve V^(A) was fitted to extend the range of A. This figure shows two remarkable features: 

1. The general trend of the constant-a curves is very similar to the trend of the experimental 
data; 

2. The a = 1 curve, which indicates the load at which the mountain-pass energy equals the 
stored energy in the prebuckled cylinder, appears to be a lower bound to the data. 

One may also consider an alternative way of rescaling energy. The cylinder is a long structure, 
and it is not clear to which extent the length of the structure is relevant for the imperfection 
sensitivity. It may be reasonable to compare the energy of the mountain pass with the stored 
energy contained in a representative section of the cylinder; the radius R provides a natural 
length scale for such a representative section. 

Similar to Figure 5.2, in Figure 5.3 we present curves of constant where (3 is the ratio of 
mountain-pass energy to stored energy in a section of length 27ri?: 

1 ^ ^(^) (52) 

^"47rV3(l-^2) R A2 ■ ^^-'^ 

Once again to obtain Figure 5.3 we fitted V{\) from Figure 4.3 to extend the range of A. 

6. Discussion and conclusions. The mathematical results and their interpretation in the 
context of a loading problem have brought a number of new and improved insights. 
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Fig. 5.2. Same data as Figure 1.1; added are two curves of constant a = 1, a = 0.1 where a is given in (5.1 ). 
Note that the load at which the mountain-pass energy equals the stored energy in the prebuckled cylinder (a = l), 
appears to he a lower hound to the data. 
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R/t 

Fig. 5.3. Experimental data and fit to Vmp with (3 = 1 and f3 = 0.1 in (5.2). Again the load at which 
the mountain-pass energy equals the stored energy in a representative portion of the prehuckled cylinder ((3 = 1) 
appears to he a lower hound to the data. 



6.1. The cylinder has doubly-localized solutions. The subcritical nature of the bifur- 
cation in Figure 1.2 strongly suggests that equihbria exist with deformation locahzed to a smah 
portion of the cyhnder length. In [20, 24, 25] such solutions are indeed calculated numerically and 
investigated analytically; these solutions are periodic around the cylinder and have exponential 
decay in the axial direction. 

The localization in the axial direction demonstrated by these solutions is consistent with results 
on simpler systems such as the laterally supported strut [21, 26]. The behaviour of the cylinder in 
the tangential direction is not as well understood. The lack of localization in the simply-supported 
fiat plate [13] suggests that the cylinder should also prefer tangentially delocalized solutions, as do 
most of the experiments. The single- and multiple-dimple solutions of this paper, however, clearly 
demonstrate that doubly- localized solutions do exist, and that some of these can be stable under 
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constrained shortening. 

6.2. The mountain pass is a single-dimple solution. The fact that the mountain-pass 
solution exists fohows essentiahy from two features, the local minimality of the unbuckled state 
and the existence of a large-deflection state of lower energy. The former is a simple consequence^ 
of the subcritical load level, but the latter is based on an essential property of the cylinder: for 
a sequence of cylinders for which R/t oo, the nondimensionalized load-carrying capacity (the 
highest load at which the unbuckled state is not only a local but also a global energy minimizer) 
decreases to zero. This property was demonstrated implicitly by Hoff et al. [15], and Lemma 3.2 
provides a simplified proof of this result and a simple sequence of functions that illustrates the 
property. 

However, the fact that the mountain-pass solution is localized, and even is the most localized 
solution that is possible — a single dimple — is interesting in its own right, and provides a comple- 
mentary view of the discussion of localization above. A different way of formulating this result is 
that "creating the first dimple is the major obstacle" ; afterwards one may increase the size of the 
dimple and add further dimples without ever returning to the same high energy level. In itself 
this interpretation points to a relationship between single dimples and imperfection sensitivity. 

6.3. Single dimples in other contexts. Interestingly, single dimples have appeared in the 
literature in a number of seemingly unrelated ways: 

• In the celebrated high-speed camera images of EBlinger [12] the first visible deformation 
is a single, well-developed dimple half-way between the ends of the cylinder. New dimples 
quickly appear next to this first dimple, and the deformation then spreads around the 
cylinder and in axial direction. It is remarkable, though, that the first visible deformation 
is a single dimple. 

• Some of the "worst" imperfections calculated by Demi and Wunderlich [9] and Wunderlich 
and Albertin [35] are in the form of a single dimple; as the load decreases, the dimple 
contracts and becomes even more concentrated. 

• Hiihne et al. [19] assert that single dimples are also realistic and stimulating imperfections 
in the sense of [34]. 

• Zhu et al. [38] base their analysis of the scaling behaviour of the experimental buckling 
load on the behaviour of a single dimple in other structural situations (the point-loaded 
cylinder and the sphere under uniform external pressure). 

Note that the single-dimple appearances above are of three different types. EBlinger's dimple 
is an experimental observation; the dimples of Wunderlich and co-workers and of Hiihne et al. are 
geometric imperfection profiles; and the dimples studied by Zhu et al. are only analogies, since 
they are solutions of different problems. 

6.4. Scale-invariance of the localized solutions. It is an interesting observation that 
the Von Karman-Donnell-equations can be rescaled to depend only on the (rescaled) load level. 
For localized solutions, for which the boundary plays no role of importance, this implies that the 
set of solutions reduces to a one-parameter family. This allows for efficient computation of the 
behaviour of such solutions but also gives interesting insight into the relationship between dimples 
in cylinders of varying geometry (see e.g. Figure 4.4). 

Naturally the scale-invariance is expected to break upon replacing the Von Karman-Donnell- 
equations by a different (probably more detailed) shell model. Nonetheless, it may reasonably 
be expected that much of the understanding of the relationship between cylinders of different 
geometries remains roughly correct. 

We certainly also expect that the large-scale geometry of the energy landscape does not depend 
on the specific model of the cylinder. Using a discrete mountain-pass algorithm to find mountain- 
pass points therefore does not depend on the Von Karman-Donnell-equations, and should give 
similar results regardless which shell model is used. 



^On a finite domain this consequence is indeed simple; on an infinite domain it appears that besides the 
third-order term also the fourth-order term in the energy has to be taken into account, as remarked in Section 3 
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6.5. Connection with sensitivity to imperfections and "Perturbation energy". 

Kroplin and co-workers [23, 11, 32] were the first to suggest an estimate of the stability of the 
unbuckled state in terms of the ratio of a "perturbation energy" (Storenergie) to the pre-buckling 
energy In early papers [23, 11] the perturbations are still fixed rather than determined, but from 
both the introduction and the final results in [32] it may be deduced that an optimization is done 
over all perturbations (although this is simultaneously contradicted on page 333 of this paper). 
Unfortunately, these papers do not provide enough details to determine exactly what the authors 
calculate. 

There is one aspect in which our method can clearly be seen to differ from these earlier 
approaches. The discrete mountain-pass algorithm takes into account global features of the energy 
landscape, and provides a global measure of the separation barrier between two states that lie 
far apart. This is different from the papers mentioned above, in which the method uses only 
local information (reflected, for instance, in the assumption that the equilibria in question lie on 
the same bifurcation branch). This difference is illustrated in Figure 6.1, where a local analysis 
might find stationary point but the mountain-pass algorithm will find the more important 
obstacle 



W3 



Wi 



W4 



W2 



Fig. 6.1. While a local algorithm to find a critical point may settle on a minor critical point such as ws^ the 
mountain-pass algorithm, by its global setup, will converge on to the essential obstacle W4. 



Appendix A. Proof of Lemma 3.2. 

Lemma 3.2 states that there exists a sequence of functions ws, 1 -periodic on , such that 

J A0^=O((52-") as 6^0, (A.l) 



[-l/2,l/2]2 [-l/2,l/2]2 

and 

[-1/2,1/2] 



for any a > 0. Here the function (j)s solves equation (2.5) with periodic boundary conditions. In 
addition, ws and (j)s satisfy (2.6) on the boundary of [—1/2,1/2]^. 

The proof consists of three parts. In the first part we construct the functions ws'^ in the second 
part we study the symmetry properties and the support of the right-hand side of (2.5); and in the 
third part we show that this sequence has the asserted scaling. 

A.l. Construction of ws- Let be given by 

f:{s) = [h dist(.,Z)<^ ^.^^ /.(0) = /aO) = 0. 

I otherwise. 

Note that / is even and that /(I) = 1/4. Define 

W5{x, y) = fs{y + x) + fs{y - x) - ^/5(2x) - ^y'^. 
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Fig. A.l. The functions f and —w; on the right the plotting area is slightly larger than one period. 

We shall drop the subscript 5 and simply write w and /. 

The function w is periodic on with period 1 in each direction. To show this, we prove that 
the first two derivatives match up on opposite sides of [—1/2, 1/2] x [—1/2, 1/2]: 

• By the symmetry of / the function w is even in both x and y. Consequently w takes the 
same values on (1/2,?/) and (—1/2,?/); the same holds for (x,±l/2). 

• For the comparison of the first derivatives we calculate 

1/2 .1/2 .1/2 .1/2 

w^^{x, y)dx= / f\y + x) + / f\y - x) dx - 2 / f\2x) dx = 0, 



1/2 J-1/2 J-1/2 J-1/2 

implying that Wx{—l/2,y) = Wx{l/2^y);hy the symmetry of w it follows that Wx{—l/2,y) =| 
Wx{l/2^y) = 0. Similarly, using the definition of fs we find that 



/1/2 .1/2 .1/2 

Wyy{x, y)dy= / f\y ^ x) dy ^ f\y -x)dy-l = {), 

-1/2 J-1/2 J-1/2 



implying that Wy{x^ ~l/2) = Wy{x^ 1/2) = 0. 
Periodicity on then follows from the remark that all second derivatives of w are periodic with 
period 1 in x and y. 

A. 2. Support, symmetry, and boundary conditions. Next we investigate the right- 
hand side of (2.5). We find 

[w, w] + Wxa: = {(/" (^ + x) + f {y - x) - 2 f " (2x)) {f " {y + x) + f'\y - x) - I) 

-{r{y^x)-r{y-x)f] 

^U"{y^x)^f"{y-x)-2f"(2x)) 
= 4r{y + x)r{y -x)- 2r{2x){f\y + x) + T (^ - x)). 

This expression has zero integral over [—1/2, 1/2]^. This follows from the periodicity of w: 

.1/2 .1/2 .1/2 .1/2 

/ / WxxWyy dxdy = / wly dxdy (A.2) 

J-1/2 J-1/2 J-1/2 J-1/2 

by partial integration. More is true, however; we analyse the support of [w, w] -\-Wxx in [—1/2, 1/2]^ 
in more detail. 

The value of f'^ is either (4(5)~^ or zero; in order to determine [w^w] + Wxx it is therefore 
sufficient to calculate the measures of the pairwise intersections of the supports of f"{y + x), 
r(^-x), and r{2x): 

• The intersection of the supports of f^\y + x) and f'\y — x) has total area 4(5^ (see 
Figure A.2). 

• The intersection of the supports of f^\y + x) and f"{2x) also has total area 4(5^ (see 
Figure A.3). 
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r 

2 



^=2 



Fig. a. 2. The areas of the black regions add up to 45^ 




Fig. a. 3. The areas of the black regions add up to 46^ 

Since the support of [w, w] -\-Wxx concentrates onto a discrete set of points, let us examine the 
behaviour at one of these points. For smah S the support forms disjoint sets in [—1/2, 1/2]^, and 
we can restrict our attention to the origin alone. 

If \s\ < 1/2, then f'§'{s) can be written as 



where 



\a\<l 
otherwise. 



Therefore, as long as \{x^y)\ < 1/4, 



fsiy + ^)fs{y -^)- 2n'{2x){n'{y + x) + n'{y - x)) 



S29 



y -\- X 



y-x 



2x 



lF{6-\x,y)), 



y -\- X \ I y — X 
9 



(A.3) 



where we introduce a new function F, which does not depend on (5, to summarize the line above. 
Note that suppF C [—2, 2]^. Note also that by (A. 2) the function F has zero integral; in addition. 
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since f'^' is even, the function 4/^'(?/ + x)f'^{y — x) — 2f'^'{2x){f'^'{y -\- x) -\- f'^{y — x)) is also even 
in X and in y. Therefore 



y)) dxdy = / yF[{x^ y)) dxdy = 0. 

This property wih be used below. 

The assertion also states that the functions w and (j) satisfy (2.6) on the boundary of [—1/2, 1/2]^. | 
We first note that w and (j) are periodic in the following sense: 

w{x±l/2,y±l/2) =w{x,y) and ± 1/2, ± 1/2) = (l){x,y). (A.4) 

For w this is a simple consequence of the functional form of w; for (j) it is a consequence of the 
uniqueness of solutions of (2.5) under periodic boundary conditions. The periodicity of w and (j) in 
the ^-direction in (2.6) then follows from a repeated application of (A.4). Similarly, the symmetry 
conditions in x in (2.6) follow from a combination of the symmetry of w and (j) around {y = 0} in 
combination with (A.4). 

A. 3. Scaling properties. We now use the information gathered above to show that the 
sequence ws has the scaling properties of (A.l). All function spaces are on [—1/2, 1/2]^. 

First, remains bounded on bounded sets as (5 — > 0; therefore w"^^ converges to a finite, 
positive value. In addition, all second derivatives of ws remain bounded in L^, so that 

\\Awsh,<C. 

The second derivative is bounded by 1/4(5, so that we can estimate 

Turning to ^5, we start by remarking that [ws^ ws] + wsxx is bounded in L-^, since 
/ \[ws,ws]^ws:,:,\ = ^ [ \F {S'^x , y)) \ = 0{1) . 

bmce W^^P ^ ^00 ah p > 1, the solution of 

satisfies 



SO that 

||^5|Ih^2,p' < C\\[ws,ws] +^^11^1 < C, 
where 1/p^l/p' = 1. Using W^y ^ ^71,1-2/^' g^^^ 

||^5||ci,l-2/p' < C \\(l)s\\y^2,p' <C. 

Writing, locally at the origin, 

Mx, y) = 05(0, 0) + V05(0, 0) • (x, y) + 0(|(x, y f^^-^'""^), 



22 



JIRI HORAK, GABRIEL J. LORD, AND MARK A. PELETIER 



we find by multiplying the equation (2.5) by and integrating, 

II^^5|Il2([-1/2 1/212) = 2 / (l)6{[ws,Ws]^WSxx} 
^[-l/4,l/4]2 

^ i[-l/4,l/4]2 

^2^^. f {x,y)F{5-\x,y))+0{5^^'-''^'^) 

J[_i/4_i/4]2 

= 0(52(i-W)), 

since the zeroth and first moments of F are zero. Since p' may be chosen arbitrary large, this 
estimate concludes the proof. 

Appendix B. Derivation of the Von Karman-Donnell equations. 

The common aim of the many elastic shell theories is to approximate three-dimensional elas- 
ticity by a reduced description in which the unknowns are functions of not three but two spatial 
variables, see for example [6]. For the Von Karman-Donnell cylinder the central approximation 
is the director Ansatz^ which states that a normal to the centre surface remains normal through 
deformation. By this Ansatz the displacement is fully characterized by the displacement vector 
{u^v^w)^ a function of the two in-plane spatial variables x and where u is the displacement 
in axial (x-), v in tangential (?/-), and w in the radial direction. Apart from some rescaling the 
function w is the same as the unknown w in the rest of this paper. 

In the formulation of Section 2 the unknowns u and v are replaced by the Airy stress function ^, 
which is derived by minimization with respect to the displacements u and v for fixed w. This 
minimization argument is well known in the context of the Von Karman plate theory and can 
be found in many textbooks. Determining the boundary conditions that the function (j) satisfies, 
however, is not straightforward (see also the discussion in [27]), and it is for this reason that we 
now describe the argument in detail. The main goal is to show that the function (j) is periodic in 
the tangential direction. 



B.l. Energy and shortening. All quantities in this appendix are dimensional. We assume 
a cylinder of thickness length L, and radius i?, and we set Vt = [0,1/] x [0,27ri?]. The stored 
energy given by [31] is 

E f t f 

" 24(1 - 1,2) ^2eJ^ ^^"'^ + " + '^)('^ll'^22 - <T?2)] 

for a linear material of Young's modulus E and Possion's ratio u. Under the assumption of plane 
stress the stress and strain tensors are related by 

l-iy^^ ^ l-iy^ V)£i2 822 + ^^11 J 

and under the small-angle approximation the strain tensor can be expressed in the displacements 
as 

_ / \wl \Uy + \Vx + \WxWy 

\\Uy + \Vx + \WxWy Vy + ^W^ " pW 

These choices for the energy and for the stress and strain tensors are very similar to those for a 
flat plate. The intrinsic curvature of the cylinder, of magnitude p = only appears in the last 
term of £22, —pw^ which expresses the fact that radial displacement creates extensional strain in 
the ^/-direction. 




CYLINDER BUCKLING 



23 



The average axial shortening is given by 

Si = [^(^' y) - ^(0' y)] dy = ^^^y^ 

and an equihbrium {u^v^w) at load level P is a stationary point of the total potential Vi = 
Ei-PSi. 

B.2. Boundary conditions. At the boundaries ^ = 0, 27tR it is natural to assume that 
and w are periodic, but at x = 0, L there is a certain amount of choice. 
The boundary conditions on w (see (2.6a)) are 

Wx = {Aw)x =0, at X = 0, L, (B.3) 

and these conditions signify a fixed angle {w^ = 0) and zero radial force {{Aw)x = 0). They may 
also be understood as symmetry boundary conditions, as in the case of a sequence of cylinders 
stacked on top of each other. For u and v we assume boundary conditions 

Uy = and cri2 = 0, atx = 0,L, (B.4) 

which signify that the ends of the cylinder are rigid in the x-direction and that there is no friction 
between the cylinder and the apparatus holding it. Note that the pair of boundary conditions 
(712 = and {Aw)x = together states that the loading apparatus exerts only axial forces on the 
cylinder. 

The boundary conditions on w are invariant under addition of a constant to i.e. under 
replacing w hy w -\- c; for stationary points we may exploit this fact. 

Lemma B.l. If (u^v^w) is a stationary point of Ei — PSi under boundary conditions (B.3- 
B.4), then 

I (722 = 0. (B.5) 



Proof Under the replacement w w -\- we have 

da _ Ep 0\ 

and therefore 

= ^{Ei - PSi) = --r^ [ [{an + a22)(^ + !)-(! + ^)(an + ^^22)] = -tp [ (122- 

□ 

B.3. Derivation of the Airy stress function cj). The energy (2.7) and the Airy stress 
function (f) are derived from the total potential Ei — PSi by minimization with respect to the 
displacements u and v for fixed w. Performing this minimization on the second term in Ei yields 
the classical plate equilibrium equations 

(7iix + cri2y = and ai2x + (722y = 0. 

Note that the derivative of with respect to u only creates boundary terms. By applying three 
times the well-known characterization of divergence-free vector fields as rotations of scalar fields 
(see e.g. [3, Th. XII. 3. 5]) we obtain the local existence of a function (j) satisfying 

(Jii = E(j)yy, (712 = -E(j)xy, and (722 = E(l)xx, (B.6) 

where we use the traditional scaling of (j) by Young's modulus. 
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B.4. Boundary conditions on cj). The existence of the function (j) is the result of a local 
differential-geometric argument, and as such gives no reason for (j) to be periodic in y. The following 
theorem shows that after a normalization transformation the function cf) can indeed be assumed 
to be periodic in and may be taken to satisfy the same boundary conditions as the function w. 

Theorem B.2. If u, v, and w are periodic in y and satisfy boundary conditions (B.3-B.4), 
then there exists a function (j) that satisfies 

(Jii - / (7ii =E(l)yy, cri2 = -E(l)xy, and cr22 = Ecl)^^, (B.7) 

is periodic in y, and satisfies boundary conditions 

0^ (A0)^ =0 atx = 0,L. (B.8) 



Remark B.l. Mechanically the normalization of (j) with j^crii means that (j) represents the 
deviation from the unbuckled in-plane stress state. 

Proof. As discussed above, there exists a function (j) satisfying (B.6); we will construct in 
stages a new function cj) which satisfies (B.7) and the boundary conditions. 
We first convert condition (B.6) into (B.7). Set 

Since the second derivatives of (j) can be expressed in terms of derivatives of ix, and the second 
and higher derivatives of (j) are automatically periodic in y. Therefore 

^P(^) = [(t^xy{x, 27tR) - (t)xy{x, 0)] = 0, 

implying that p is actually independent of x. (A mechanical argument provides the same result: 
Etp{x) = tj^Q^^ aii{x^y)dy is the total force applied at a virtual cut at level x, and mechanical 
equilibrium implies that this force is independent of x). Therefore 



ft\p = 27rR f pdx = [ <l)yy = [ ^^11, 
Jo Jn ^ Jn 



so that the new function 



~ p 
(j){x,y) := (l){x,y) - -y^ 

satisfies (B.7). Note that this implies 

lyy = 0. (B.9) 

We now turn to the periodicity in the ^/-direction. It remains to show that (^, (px^ and (j)y are 
the same at = and y = 27tR. Again the periodicity of the second derivatives implies that 

— [(t){x, 2iiR) - (f){x, 0)] = (pxxix, 27iR) - (pxxix, 0) = 0, 

so that ^(x, 27rR) — ^(x, 0) = ax -\- b for some a, 6 G M. Defining 

7 / \ 7 / \ ^y / / \ P 2 ^y 
0(x, y) := 0(x, y) - —— = 0(x, y) - -y 



27iR 2^ 27iR 



CYLINDER BUCKLING 



25 



the function still satisfies (B.7), and 

(j){x^ 27tR) — (j){x^ 0) = ax. 

Finally, we find that 

a = (^,(0,2^i?)-(^,(0,0) = ^:,y{0,y)dy = - ai2{0,y) dy ^ = ^ 0, 

and therefore that (j){x^27TR) — ^(x,0) = for all x. The same follows for (j)x{x^27TR) — (j)x{x^O) 
by differentiation. 

To show that (j)y also matches, 

^ [^y{x, 2tiR) - ^y{x, 0)] = [(t)xy{x, 2iiR) - (t)xy{x, 0)] = 0, 

and therefore (j)y{x^27rR) — (j)y{x^O) is constant in x; by (B.9) this constant is zero. This proves 
that (/) satisfies (B.7) and is periodic in y. 

We finally discuss the boundary conditions at x = 0, L, and we follow the line of reasoning 
of [27]. By (B.4) and (B.l) £12 = at x = 0, L, so that by (B.3) and (B.4) 



d_ 

dy 

Therefore 



xy — ^(2Sl2 — Uy — WxWy) =0 at X = 0, L. 



£22x = Vxy + U)yWxy - pw^ = at X = 0, L. 



Using Ee22 = cr22 — ^cni, we then find 

1 d 

E d~x 



^(f^xyy — 4^xxx — ^4^xyy — ——{(^22 — ^(^ll) — ^22x — 0, 



and by adding (1 + i^)4>xyy = ~(1/^)(1 + T^)^i2y = it follows that 

(A(^)^ = 4>xxx + i>xyy =0 at X = 0, L, 

which proves one part of (B.8). 

From (j)xy = —cri2/E = we find that 

^x(0, y) = Co and 4>x{L, y) = cl for ah y G [0, 27tR]. 

Writing 

27rR{cL - Co) = / 4>xx = ^ I cr22 ^ 
we find that cl = cq. Now the function 

^(x, y) := (^(x, y) - cqx = y) - |?/^ - - cqx 
satisfies (B.7) and (B.8) and is periodic in y. This concludes the proof. □ 

Remark B.2. It is instructive to note that the periodicity of (/) is a result of the specific choice 
of boundary conditions, and will in fact not hold if different boundary conditions are taken. For 
instance, if a tangential shear stress r is applied at the cyhnder ends {i.e. the cylinder is loaded 
under torsion), then the coefficient a in the derivation above will not vanish, and (l)x will not be 
periodic in y. 
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B.5. Putting it all together. By an elementary but lengthy calculation we find that ^, as 
provided by Theorem B.2, satisfies the equation 

A'^(l)^ pwx:^^[w,w]=0 inn, (B.IO) 

and that the second term in Ei can be written as 

A(/)2 -2(l + ^)[0,( 



tE r 



By the boundary conditions given by Theorem B.2 the second term vanishes, and the total stored 
energy functional can therefore be written as 



Note that this energy is a function of w alone; the function (j) in this definition is assumed to be 
given by (B.IO), with the boundary conditions of Theorem B.2. Similarly, we rewrite the average 
shortening as 

S2{w) := Si{u) = --^ / 



27iR 



(B.2) 



(B.l) 1 /■ r 1^1 

(B.7) 1 f r. , n 1 

= -2^yj^--^^-]+4^., 



1 



AttR 

A stationary point w of E2 — PS2 satisfies the Euler equation 

f2 p 



12(1 2i^REt 
where again (\) is related to w by (B.IO). With the nondimensionalization 

w = 4:7r^Rw, (j) = Wtt^R^ (j), x^2tiRx, y^2TiRy, 

we then obtain equations (2.4) and (2.5), and the dimensional energy E2 and average shortening 
5*2 above can be expressed in these variables as 

E2 = ^^^^^'^2^ J + 327rhER^ J a/, ^'2 = ^tt^'R J wl (B.ll) 
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